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The time evolution of structure factors in the disordering process of an initially phase separated lat- 
tice depends crucially on the microscopic mechanism which drives the disordering, such as Kawasaki 
dynamics or vacancy mediated disordering. Monte Carlo (MC) simulations show unexpected "dips" 
in the structure factors, which mean-field theory completely fails to capture. The disordering via 
vacancies is slower by a surprisingly large constant factor compared to Kawasaki dynamics. To 
understand the dips in odd structure factors (i.e., those reflecting the initial order), a phenomeno- 
logical model is introduced, and an analytical solution of Kawasaki dynamics is derived, in excellent 
agreement with simulations. The presence of dips in the even structure factors for vacancy mediated, 
and their absence for Kawasaki dynamics, marks a significant, but not yet understood, difference of 
the two dynamics. 



After a temperature "downquench" into a coexistence 
region, binary systems such as magnets or alloys develop 
long-range order in a universal fashion jl]. Determined 
by symmetries and conservation laws, the dynamic coars- 
ening follows one of a small number of scenarios, char- 
acterized by dynamic scaling and universal exponents. 
Considerable work has been devoted to the theoretical 
foundations of this behavior By contrast, far less 
attention has focused on the opposite phenomenon, i.e. 
disordering following an upquench in temperature, even 
though it is hardly less important in materials science: 
It determines, e.g., the time evolution of an interface be- 
tween two materials, or the waiting times which must 
elapse before a heated material has mixed sufficiently. 

A key ingredient in the physics of these processes is 
the kinetic mechanism by which two particles exchange. 
In many materials, it is dominated by the diffusion of va- 
cancies or defects 1^, rather than direct atom-atom ex- 
changes. Vacancy concentrations are typically very small 
10~^) 1^. For this type of dynamics, it has recently 
been shown |^ that a simple mean-field theory suffices 
to obtain a quantitative description of disordering charac- 
teristics which are controlled by short-range correlations. 

In this letter, we will show that this mean-field the- 
ory fails dramatically, even at the qualitative level, to 
predict the dynamics of observables which include cor- 
relations over larger distances in the system. A typical 
example is the time-dependent equal-time structure fac- 
tor. We present MC simulation data and a microscopic 
theory, illustrating the crucial role played by the details 
of the whole vacancy path. Comparing data for vacancy- 
mediated and direct exchange dynamics, we obtain a key 
result of our study, namely, that experimental data for 
the equal-time structure factor can identify the micro- 
scopic kinetic mechanism by which a material evolves. 
We conclude with a list of open questions. 

We consider a c?-dimensional square lattice of size L"^ 



{L even, for simplicity) with periodic boundary condi- 
tions. Each site ? is occupied by one of two types of 
particles, which we label as spins: sj; — ±1. A single 
spin is tagged, playing the role of a vacancy or defect. 
Only this "vacancy spin" is allowed to move, all others 
being passive. Turning to the microscopic dynamics, we 
define two versions: (i) The vacancy performs a Brown- 
ian random walk on the lattice, exchanging with nearest- 
neighbor spins and thus scrambling the spin configura- 
tion, or (a) the tagged spin performs one single step, and 
then passes the tag to another, randomly selected spin. 
This sequence is then repeated. Note that both versions 
strictly conserve particle number (or total magnetization, 
sr). The first is known as vacancy mediated dynam- 
ics (VMD) l3|j], the second is precisely Kawasaki dy- 
namics (KD) In the language of Ising models, they 
correspond to T = cxd dynamics, since energy costs are 
ignored. Yet, we consider a perfectly phase-segregated 
initial configuration (up- and down-spins separated by 
two planar domain walls perpendicular to, say, the x- 
axis), corresponding to T = 0. At time t — 0, the system 
experiences a temperature upquench, from zero to infin- 
ity, and we monitor the disordering process as a func- 
tion of time. For t ^ oo, the system clearly reaches the 
infinite temperature equilibrium state, characterized by 
completely random configurations. 

We monitor the time-dependent structure factors, 

Sdt) ^ L-'^^ 



X;,-sj-(t) exp(-ifc- j) 



Here, (•) de- 



notes the time-dependent average over a large number 
(1,000 for our data) of independent runs, starting from 
the same T — Q configuration differing only in the ran- 
dom initial position of the tag. MC time is incremented 
by 1 for each exchange. The most interesting structure 
factors are those whose wave vector k is perpendicular to 
the initial phase boundary, k = (27rn/L, 0, . . . , 0), with 
n integer. For odd n <^ L, they are sensitive to the 
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initial order so that Sj:{0) is of the order L"^, while for 
even n, ^(O) = (since L is even). For brevity, these 
two classes will be referred to as odd and even k. Clearly, 
linif^oo ^/^{t) = 1 for all fc, due to the random final state. 




FIG. 1. Time evolution of the first eight structure factors 
with k = (27m/i,0), n = 1,3,5,7,8,6,4,2 from top to 
bottom, on a d = 2 lattice of size L = 160 with (a) VMD 
and (b) KD. The thick solid straight lines in the insets, which 
show the magnified crossover region, have slope 1/2. The two 
dashed vertical lines window the intermediate regime. 

For reference, we review a few key results from The 
disordering process displays three regimes, separated by 
two characteristic times, tE — L^, and ii, ~ L'^'^^. For 
t < tE only a fraction of sites has been visited by the 
vacancy, and the equilibrium state is reached for t > i^- 
Between these well-separated bounds lies the "interme- 
diate regime," in which the vacancy visits each site many 
times without destroying much of the initial order. Here, 
the disorder parameter, being the number of broken near- 
est neighbor bonds (NN), displays dynamic scaling with 
universal exponents before reaching its equilibrium value. 

Turning to structure factors, in Fig. |l| we compare MC 
data for a 160 x 160 lattice, with (a) VMD and (b) KD. 
First, note that three regimes emerge again. The final 
regime is easily identified by all S'^'s reaching unity. In 
the intermediate regime, the odd S^'s clearly display the 
decay of the initial order which remained largely unaf- 
fected during the earliest regime. Second, we emphasize 
a rather peculiar feature, namely, that several structure 
factors develop a minimum, manifested as a small "dip" 



before they finally tend to their equilibrium value. For 
VMD this dip appears in all structure factors except the 
two lowest ones, i.e., n = 1, 2 (strictly speaking, in n = 4 
it emerges as a "shoulder"). In contrast, for KD it exists 
only in the odd Sj^'s. At first sight, these dips might sug- 
gest a temporary rebound towards order. We will show 
below, however, that their origin is different. 
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FIG. 2. Short time behavior of the first 4 even structure 
factors on a d = 2, L = 500 lattice with VMD, averaged over 
1,000 runs. The thick solid line has slope 1, the vertical line 
indicates the beginning of the intermediate region. 

Third, turning to "short" times, just after the onset of 
the intermediate regime, we observe a clear difference in 
the even structure factors for the two types of dynam- 
ics, displayed more visibly in the data for a much larger 
(i — 500) lattice (Fig. The short-time behavior for 
a VMD system appears to be ^ with a short- 
time exponent 7v ~ 1, whereas the KD system follows 
Sj: ^ with 7k > 1. Below, we will see that 7k — 3/2 
exactly, in the i — > oo limit. Thus, the presence or ab- 
sence of dips and the short-time behavior offer a simple 
signature to distinguish VMD from KD. Finally, Fig. |l| 
indicates that disordering by KD is faster than VMD by 
a factor of about 4.5. This factor is roughly independent 
of the lattice size but decreases with d 0| . We emphasize 
that, in all figures, the error bars are of the order of the 
line thickness and therefore not shown. 

We now turn to the theoretical descriptions of our 
model. Our first approach, an intuitive model of two 
"fluids," is based on the coarse-grained mean-field the- 
ory of and accounts for the dips in the odd structure 
factors. The second method involves an exact equation 
of motion for the microscopic spins sj{t)^ which can be 
solved for the KD case and leads to 7k — 3/2. 

A "two-fluid" model. To motivate this approach, we 
briefly recall the key ingredients of the mean-field the- 
ory 1^. Coarse-graining the microscopic dynamics in 
space and time, we identify two slow variables, namely, 
the local vacancy and magnetization densities, Lp{r, t) 
and il){r,t). In the intermediate regime, the vacancy 
distribution has already reached its equilibrium value, 
tfo = L^'^, so that the only non-trivial dynamics is 



2 



embodied in ■0(r, t), which obeys a simple diffusion 
equation: dtip = DV'^tp. Here, the diffusion coef- 
ficient D (X ipo reflects the presence of a single va- 
cancy. The solution, with initial condition ^(r, 0) = 
sgn(x), |a;| < i/2, and periodic boundary conditions, 
is i^{x, t) = Y.k llLe*'=^e-^'='*5„ odd, where k = 2Tm/L 
and 5n odd = 1 for n = ±1, ±3, . . . and otherwise. 

We propose to describe the local spin density by the 
incoherent sum of two components ("fluids"), represent- 
ing the ordered and disordered fractions. The ordered 
component, starting from being 100% phase segregated 
and ending at zero, can be taken as ip^Xjt). The remain- 
ing component , with density 1 — | -0 (a; , i ) | , is fluctuating so 
that it will be represented by (1 — \'ip{x,t)\)ri{x,t), where 
r]{x,t) is a delta-correlated noise with zero mean. Sum- 
marizing, we write the spin density s{x,t) as 



s{x,t) ^ i^{x,t) + {1 - \i^{x,t)\)f]{x,t), 



(1) 



Taking the Fourier transform of s(a;, t) and averaging over 
rj to obtain Sk{t), we find the final result 



Skit) = L"" 



^71 O 



dd 



E 

n' odd 



(2) 



where r = L'^/{4tt^D) plays the role of ti,. The first 
term captures the decay of the initial ordered fraction. 
Note that only the odd fc's are present. By contrast, the 
second term is k-independent and reflects the buildup of 
the disordered component. 

The origin of the dips is now apparent: The first term 
starts with a macroscopic amplitude, but decays on the 
time scale Tjlv?. Meanwhile, due to the sum over all 
n', the second term rises to unity with r, regardless of 
n. Therefore, Sk{t), for odd fc, does not decay to unity 
monotonically. Further, these dips occur at earlier times 
for larger fc's, a feature which agrees with the data, at 
least qualitatively (Fig. |l]) . We believe that the contrast- 
ing terms may be interpreted physically, the first being 
associated with the relatively fast deterioration of the 
initial phase boundary while the second represents the 
slower buildup of the disordered component. 

Eq. (H) leads us to another prediction, namely, how 
the final state is approached. By Poisson resummation 
H, it can be shown that the second term increases as 
i for times t < t. This effect is most easily seen in 
the even S^^s where the ordered contribution is absent. 
The simulation data show excellent agreement with this 
behavior for both VMD and KD (Fig. g). 

Not surprisingly, this simple approximation has its 
shortcomings, such as J2k^k{t) L''-. Also, according 
to Eq. (||), all even structure factors share the same be- 
havior, which is clearly not borne out by the simulations. 
Moreover, Eq. (|^) gives no indication of an early-time 
crossover from a t''' to a i^/^ behavior. To improve our 



understanding, we now turn to an exact formulation of 
the dynamics. 

Exact lattice calculation. In the following, we analyze 
the microscopic motion of the vacancy on the lattice. The 
resulting equation of motion for the structure factor is, 
of course, exact. Remarkably, it can be solved in the 
case of KD, providing us with detailed information about 
different temporal regimes and finite size effects. Here, 
we summarize the key features of this approach, while 
deferring all details to a later publication |Q. 

Clearly, the dynamics at each time step (for both 
VMD and KD) depends on the exchange of only two 
spins, so that the Liouvillian in the master equation 
dtP{{sj}-t) = Y.is'^C. {{st};{s',})P{{sL}-t) consists of 

a sum over terms of the form (5(s^; sj)J(sr; Sj). As a re- 
sult, the equation for G{x]t), the two-point function, is 
linear and involves no higher correlations. Specifically, in 
the KD case, one obtains (for x ^ ) 



—dtG{x;t) 



AjG(f;i) 



X^%[G(6;t)-G(0;i) 

{b} 



(3) 



where dt, Aj, and {b} stand for discrete time differences, 
the lattice Laplacian and the set of 2d lattice vectors, 
respectively. Of course, we may replace G(0; t) by unity. 
Taking the Fourier transform of Eq. (^, we may write 
an equation for the structure factor in the form: 



(4) 



An alternative approach starts with an equation of mo- 



tion 



dtsr{t) 



V is the position of the tag at time i, and 6 is a lattice 
vector which specifies the direction of the exchange. In 
Fourier space, sjXt ^ 1) = Mj,j,,{v,h)s-^,(t), where 

^hk' ^ hk' - c-''(*^-'^')*^(l - e-*^^)(l - e^^'^)/U\ Since 
the structure factor is defined by S-^{t) — {Sf:{t)s_j:{t)), 
where the average is over all possible moves over the en- 
tire time period [0, t] , we obtain, formally, 

{p^},m 

•••A^PiPo^'%fo)%o(0)%o(0) • (5) 

For VMD, the allowed moves are correlated leading to a 
non-trivial average. For KD, however, the position of the 
vacancy is randomized before each exchange, so that the 
above average factorizes into t factors. Cast in the form 
of an evolution equation, the result is precisely Eq. (jj). 



with Xj:p 



Mki'M^k^k' 



hk' 



With either approach, thanks to the simplicity of X, 
the eigenvectors and eigenvalues can be found exactly 
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in the L ^ oo limit. Deferring details to ||7|, we quote 
only the result for d = 2. With the reminder that k — 
(27rn/L, 0) and the initial value ^(O) — {2L/T:nf 6n odd, 
we obtain 




FIG. 3. Comparison of Monte Carlo data (circles) for KD 
(same data as in Fig. |l] (b)) and the exact calculation (solid 
lines) for the same structure factors as in Fig. |l|. The dashed 
line has slope 1, the solid lines have slope 3/2 for small times. 



As in the "two-fluid" model, the two terms carry the 
same interpretations, confirming a key mean-field result: 
the associated crossover time scaling with L**. In ad- 
dition, comparing the first terms in Eqs. (0,^, we may 
identify the diffusion constant D as 1/ i^Lrj. This value 
may be interpreted as the rate of spin-exchange being 
equal to the hole density (i^^) and the two indepen- 
dent directions of its moves. The main improvement is 
the ^-dependence in the disordered component. In more 
detail, Eq. (||) correctly accounts for the presence (ab- 
sence) of dips in the odd structure factors for n > 1 
(n = 1). Turning to early behavior within the interme- 
diate regime, we may extract remarkably detailed infor- 
mation. For short times {t ~ L^), the behavior of the 
sum in Eq. (^ is controlled by the large I contributions, 
whence P — n? ~ P. Using Poisson resummation, we 
find S^[t) oc r^t^l"^ . Both the t- and the n-dependence 
are borne out by the data. For larger times, r^tjL'^ ~ 1, 
a third power law emerges which is particularly notice- 
able for large n. Here, the sum is dominated by smaller 

so that — P ~ ri^. Resummation now yields t^^^, 
with an n-independent amplitude. For large n, this power 
law sets in at earlier times, so that structure factors with 
larger n merge before being joined by Sg's with smaller 
n. Fig. H shows a direct comparison of MC data with 
Eq. (|^). Excellent agreement is observed over many or- 
ders of magnitude without any fit parameters. 



Conclusions. Using MC simulations and analytic ar- 
guments, we analyzed two types of microscopic kinetics 
(vacancy-mediated and Kawasaki dynamics). The struc- 
ture factors display several remarkable features, which 
characterize different temporal regimes. Just before the 
onset of equilibration, the competition of surface- vs. 
bulk-dominated time scales generates unexpected min- 
ima in several structure factors. These are particularly 
prevalent in the "odd" structure factors that are sensi- 
tive to the initial order. In contrast, the "even" structure 
factors grow via a sequence of different power laws before 
saturating. This sequence includes a t^^^ regime for KD 
and a t regime for VMD. Moreover, VMD generates dips 
in even structure factors which are absent in KD. All of 
these features carry over to larger system sizes and d = 3; 
in fact, the range of the power law behaviors increases 
significantly with system size. Thus, the structure fac- 
tors carry a clear signature of the underlying microscopic 
dynamics. We believe that this could lead to a simple ex- 
perimental identification of the dominant kinetic mecha- 
nism in a material, e.g., via small-angle X-ray or neutron 
scattering. 

Two complementary theoretical approaches - a phe- 
nomenological coarse-grained two-fluid model and an ex- 
act lattice calculation - provide quantitative insight into 
the physical origin of these features. Several questions 
remain open. First, our exact lattice analysis can be car- 
ried out only for KD, where the dynamics of the spins 
is Markovian. For VMD, in contrast, though the va- 
cancy path is random, it induces non-trivial correlations 
in the spins. Work is in progress to extend our analysis 
to this case [Q| . Presumably, these correlations are also at 
the root of the observation that VMD equilibrates more 
slowly than KD, by a factor of about 4.5 in d = 2. It is not 
clear at this stage, however, whether this this effect can 
be captured by an appropriate generalization of correla- 
tion factors jl^. Finally, the effect of several vacancies 
or tags should be investigated, as well as upquenches to 
temperatures other than infinity. These cases have so far 
only been studied at the mean- field level |]. 
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